RATE OF CONVERGENCE FOR DISCRETE 
APPROXIMATION OF OPTION PRICES 
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Abstract. In this article, we provide a method how one can obtain 
exact formula for the error in discrete approximation of option prices in 
continuous model. In particular, we show that it is sufficient to study 
the error for call options and digital options in order to give an explicit 
formula for the error of the price approximation of a general option. 
As an example, we study the convergence of prices of general class of 
binomial models to the Black-Scholes prices. Barrier options in CRR 
model are also considered. 
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1. Introduction 

In this article, we study discrete approximations of a continuous model 
and provide a method to obtain exact formula for the pricing error for general 
class of options. Moreover, we show that the rate of convergence is the same 
as for digital options. In addition, we show that if the payoff function is 
continuous and the second derivative exists almost everywhere, then the 
rate of convergence is the same as for call options. This also holds for 
convex payoff functions, for which the second derivative exists in the sense 
of distribution. We illustrate our result by approximating Black-Scholes 
prices with generalised class of binomial models and find the exact formula 
for the error. We also consider barrier options in CRR model. 
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In our method, we derive the exact coefficients for the error expansion 
in terms of error expansion for digital and call options. Our main tool is 
the pricing formulas derived in |15j . in which the author provides a pric- 
ing formula for general class of models and payoff functions. Moreover, 
the pricing formulas are valid for all kinds of underlying assets, including 
path-dependent options. Hence our method can be applied to obtain error 
expansion for path-dependent option as well. 

In a binomial model with n time periods, the stock price S either rises to 
Su or falls to Sd at each time period. With different choices of u and d one 
gets different kind of binomial models. 

It is well known that the binomial option prices converges to the Black- 
Scholes option prices as the number n of time periods increase to infinity. 
The first proofs to this result were given by Cox et al. [J| and Rendleman 
and Bartter [13]. A more general proof was afterwards provided by Hsia [Tj- 

Later the focus has been on the rate of convergence. Heston and Zhou [6j 
worked with Cox, Ross, and Rubinstein (CRR) model, where u = e a ^~^ and 
d = 1/u with a corresponding to the volatility of the stock in Black-Scholes 
model and At = ^. They showed that the error is of order O f° r 
general class of options. However, for European call options the convergence 
is faster. Leisen and Reimer [TTJ proved that for the CRR model, the Jarrow 
and Rudd [S] model and the Tian [T3] model, the error for European call 
options is of order O (n ). However, Leisen and Reimer did not give an 
explicit formula for the error. 

First two results on the exact formula for the error are due to Francine 
and Marc Diener [5] and Walsh |16| . Francine and Marc Diener worked with 
the CRR model and prove the following result in our notation: 

In the rz-period CRR model, if Sq = 1, the binomial price ^By( n )(-^0 °^ 
European call option with strike K and maturity T = 1 satisfies 

vSrUK) = V§ s( K) + A L ^-^K-P tl 
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where VS S (K) is the corresponding Black-Scholes price, 



A n = 1 — 2frac 



BS 

^)+nlogd 



A = -a 2 (6 + d\ + d\) + 4(d? - d|)r - 12r 2 



ogu— logd 

and frac(x) denotes the fractional part of x. 

Walsh chooses u = g^v^+r-At anc j ^ _ e -a\AAt+rAt an( j considers a gen- 
eral class of options in a binomial model with even number of periods. In 
particular, Walsh considers payoff functions which are piecewise C 2 , regular, 
and the function and its derivatives are polynomially bounded. Walsh also 
gives an explicit formula for the error. For call option he obtains 

BT{n)\ I BS\ I 24(J ^f n V«V^y 

where A is constant and 

A„ = l-2/r4' 0g ,^ + | "'° g ''' 
log u — log a 

The A n appearing in these formulas oscillates between —1 and 1, and 
thus the binomial prices usually oscillate around Black-Scholes prices. As a 
result, it can happen that doubling the number of time steps leads to larger 
error. The smooth convergence, in which case one can improve the rate by 
standard extrapolation techniques, was studied by Chang and Palmer [3]. In 
particular, they generalised the results of Diener and Diener [5] and Walsh 
|16j by considering a generalised class of binomial models with additional 
parameter A n appearing in u and d. As a result, they obtained an explicit 
formula for the error for digital and call options in their generalised class 
of binomial models, and with suitable choices of A n they obtained smooth 
convergence. The result of Chang and Palmer is one of the main tools of 
our analysis, and we introduce their result in section [H 

The rate of convergence can also be improved by constructing the binomial 
model in suitable way. The problem was studied first by Joshi [TU] for odd 
number of time periods and later by Xiao [17] for even number of time 
periods. However, in such cases one needs to consider the underlying option 
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for which one wants to improve the convergence and the improved rate does 
not hold for general options. 

In this paper, we consider generalised class of binomial models introduced 
by Chang and Palmer and derive an explicit formula for the error for general 
class of options. We consider payoff functions which are piecewise C 1 and 
polynomially bounded, and show that if the initial stock price is So, then 
the binomial price v£ T( with n time periods satisfies 



where V^ s is the Black-Scholes price of the option. We also give explicit 
formulas for coefficients Z Ut i and Z n p- Our result covers wider class of 
options than the result of Walsh. Moreover, our formula seems to be a 
simpler one and has a simpler proof. Our result is also valid if n is odd 
number. We also give additional assumptions under which we obtain rate 
of order O (n^ 1 ) . 

The fact that binomial prices converges to Black-Scholes prices holds also 
for path-dependent options (e.g. see PQ). For barrier options, the conver- 
gence was studied by Lin and Palmer [12] , who derived explicit error formulas 
for barrier call options in CRR model and showed that the rate of conver- 
gence is of order O (-^) . With our method, we obtain error expansion for 
general payoff functions. 

Our method can be used to obtain error formulas for any discrete approx- 
imation of continuous time model. For example, one can obtain error for- 
mulas in trinomial model approximations for Black-Scholes prices. Another 
example is a discrete approximation of exponential Levy model. Cawston 
and Vostrikova [2] considered an exponential Levy model and derived suf- 
ficient and necessary conditions under which the prices of certain discrete 
time Levy model converges to the actual prices. However, they did not con- 
sider the rate of convergence. Third example is the article by Xu et al. [9], 
who considered the numerical approximation for a European option pricing 
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model with jump-diffusion. In these cases, it is sufficient to derive the error 
for call and digital options and by following our method one can obtain the 
error for general options. 

The rest of the paper is organized as follows. In section [2] we state our 
main results with discussions. Section 0] contains all auxiliary facts which 
we need to prove our main results. Section [5] is devoted to actual proofs. 

2. Results 

We use the following notation: a as the volatility, r as the continuously 
compounded interest rate and T as the maturity. We denote by Vg T ,da) 
the price of the call option with strike a and maturity T in the n-period 
binomial model and Vg S (a) the Black-Scholes price of the call option with 
strike a and maturity T, given by: 

V£ s {a) = 5 *(di) - ae- rT <S>(d 2 ), 

where $(■) is the standard normal distribution function. Similarly, the prices 
for digital options l x > a with strike a are denoted by Vg T , n Ja) and Vg S {a). 
Prices for option / (St) are denoted by V^ T ,-, and V^ s . We also need several 
short notations for different functions of n and a. The list and definitions are 
given in appendices. We consider the following generalised class of binomial 
models introduced by Chang and Palmer [3J. 

Definition 2.1 (The class of binomial models). Let At = ^ and X n an ar- 
bitrary bounded function of n. We consider n-period binomial model, where 

u _ e aVM+\ n cr 2 At ^ _ e -aVAi+X„a 2 At 

with initial stock price Sq. 

With different choices of A n one gets different binomial models. For ex- 
ample, the choice X n = gives the Cox, Ross and Rubinstein model, and 
the choice X n = r/a 2 — 1/2 gives the Jarrow and Rudd model [SJ. 
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We consider the following class of payoff functions. 



Definition 2.2. For a function f : IR+ — > IR, we denote f G IT if f is con- 
tinuously differentiable except on a set of points (which is at most countable) 
< si < S2 < ■ ■ ■ < sn (and possibly on so = 0) in which f and f have 
jump- discontinuities. Moreover, f is polynomially bounded i.e. there exists 
positive constants c\,c-i and p such that \f'(a)\ < c\ + c?,a v for every a > 0. 



In what follows, the jump from left at zero and left-limits at zero are set 
to zero i.e. A_/(0) = /(0-) = /'(0-) = 0. 
Now we can state our first main theorem: 



Theorem 2.1. Let f G II. Then for the n-period generalised binomial 
model, the price of a European option /(St) satisfies 



(2.1) 



-rT -=1 

e e 2 



1 /*°° 

Vlr^^ls^i f(a) St. 



A n da 



1 



-rT -=2- 

e e 2 



f-/ /'(« , 
n Jq V27r 



Br, 



doAl 



da 



+ Yl &-f(sk)Jn(s k )+ Yl A + /( Sjfc )J n ( Sfe ) 
Ni(s k ) N!(s k ) 

+ [^-/(s k ) + A+/(Sk)]Jn(Sk) 

N 2 (s k ) 

+ o(n- 1 ), 



where the functions A n , B n , J n and J n are defined in Appendix A, N±(sk) 
is the set of jumps of f on the lattice points and AT 2 (sfc) is the set of jumps 
of f not on the lattice points. Moreover, the integrals are bounded in n. 



Theorem 2.2. Let / G II and assume that /'(a) is absolutely continuous on 
every interval (s&, Sfc+i). Moreover, assume that f" is polynomially bounded. 
Then for the n-period generalised binomial model, the price of a European 
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option /(St) satisfies 



= V BS + ~ f"(a)I n (a)da 



1 



BT(n) 



+ Y &-f(sk)Jn{Sk)+ Y A+/( Sfe )J ri ( Sfc ) 



Ni(s k ) JVi(s fc ) 



(2.2) 



+ Y, [A-f(s k ) + A + f(s k )]J n (s k ) 



N2{s k ) 



+ ~Y l (f , (s k +)-f , (8 k -))I n (8 k ) 



k=0 



+ o{n x ), 



where the functions I n , J n and J n are defined in Appendix A, Ni(s k ) is the 
set of jumps of f on the lattice points and N2(s k ) is the set of jumps of f 
not on the lattice points. Moreover, the integral is bounded in n. 

Remark 2.1. For general convex functions f , the second derivative f" exists 
in the sense of distribution as a positive radon measure fi. In this case we 
obtain rate of order 0(n^ 1 ) under additional hypothesis on the measure fi. 
This can be obtained by approximating convex function with a sequence of 
C 2 functions. 

Our theorems show that for a wide class of options, the rate of convergence 
is of order 0(-4=). Moreover, the derivative /' of the payoff function is 
usually absolutely continuous and hence the second derivative exists almost 
everywhere. In this case, the rate is of order 0(n~ l ) and if the function has 
discontinuities, then the rate is of order 0(l/ v / n). 

Walsh also gave an explicit formula for the error and proved that under his 
assumptions, the rate of convergence is of order 0(1/ y/n) if the discontinuity 
is not on a lattice point and the rate of convergence is of order 0(n~ r ) 
if all the discontinuities are on lattice points. Compared to the formula 
of Walsh, our result covers a wider class of payoff functions and a wide 
class of binomial models. Note also that we obtain, as Walsh did under his 
assumptions, that if /' is absolutely continuous on every interval (s k , s k +\), 
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all the discontinuities lie on lattice points and the function is regular, then 
the rate of convergence is of order 0(n _1 ) for any binomial model within our 
class. Indeed, if all the discontinuities lie on lattice points, then A n (sk) = 
1 for all jump points s^ and by regularity, A_/(sfc) = A + /(s^). Thus 
the coefficients of the term ly/n cancel and we obtain convergence of order 



With same procedure, we can obtain explicit error expansion for Barrier 
type options if the error formulas are known for barrier digital options and 
barrier call options. Using result of Viitasaari |15j and result of Lin and 
Palmer [12] (see Corollary 14.31 and Theorem 14. 4p we can obtain explicit error 
formulas for Barrier options in CRR model for continuous payoff functions. 



Theorem 2.3. Let f £ H and assume that f is continuous and f'(a) is 
absolutely continuous on every interval (sfc,Sfc-fi). Moreover, assume that 
H 7^ s/, Vfc € IN. Then for the n-period CRR model with initial stock 
price is So, the asymptotic formulas for European options with barrier H 
and maturity T are as follows: 

(1) for down-and-out (do) and down-and-in call (di) options with So > 



Q{n~ l ). 



H: 




1 

+ - 





k=0 




k=0 



+ o{n~ l ) 
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f£x„)(«) = ^if (a) - ^ ^ f"(a)F nA (a)da 

+ - / /"(a)F„, 3 (a)da 
n Jo 

1 N 

-7=E(/ , ( s *+)-/ , ( s *-)) F ».i(**) 

v n fc=0 

1 N 

+ - ^2(f(s k +) - f'(s k -))F n , 3 (s k ) 



n 

k=0 



oin^ 1 ) 



(2) for up-and-in (ui) and up-and-out options (uo) with So < H: 

1 f°° 

V mn)( a ) = V ls\ a ) + ^J o f"(a)F nA (a)da 

1 r°° 

+ - / /"(a)F n)5 (o)do 
n Jo 

1 N 

+ -= £(/'(* fc +) - f'(s k -))F nA (s k ) 

1 N 

k=0 
+ o(n- 1 ), 

= V^°(a) - -= / f"(a)F nA (a)da 
v n Jo 

1 f'°° 

+ - / f"(a)F nJ (a)da 
n Jo 

1 N 
l w 

+ -£(/W)-/ , (**-))* , n,7(*k) 



y ir(n)( a ) 



n 

k=0 



+ o(n x ). 

where the functions F n ^ are given in Appendix B . 



Remark 2.2. On the theorem we assumed H 7^ s k for every k € IN. 7/ this 
is not the case however, the only effect is that the formulas become slightly 
more complicated as we have to consider whether H = s k is a terminal stock 
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price or not. Indeed, the error formulas for down-and-in and down-and-out 



Remark 2.3. As for stock options, we could obtain error formulas for bar- 
rier options with convex payoff as well. 



In this article, we have derived the rate of convergence and exact formula 
for the convergence of the binomial tree scheme. We have covered a general 
class of binomial models, and our result remain valid for wider class of payoff 
functions than those of Walsh. Our formula is also quite simple as the error 
can be expressed only in terms of Lebesgue integral for continuous payoff 
functions. For discontinuous payoff functions, we only have to add terms 
that correspond the jumps. We also obtain error formulas for barrier type 
options in CRR model approximation for Black-Scholes prices. 

A significant benefit of our results is that the same method can be used to 
obtain rate for any discrete time approximation of continuous time model. 
In particular, in any such case it is sufficient to compute the error for call 
options and digital options, and the error for general options follows from our 
procedure. Note also that the same method can be used to obtain error for 
path-dependent option if the error for corresponding digital and call options 
are computed. 

As a drawback, we lose rate in the remaining term. In particular, usu- 
ally the error expansion for digital or call options are obtained by Taylor 



mial models. Following our method, the remaining term is not O ( — 7= I 
but rather o(n _1 ). However, this can be overcome by computing the ex- 
act coefficients for higher order terms in error expansions of digital and call 
options. 



options are different in these cases (see Theorem \4-4\ )- 



3. Conclusions 



expansion and the remaining term 
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4. Auxiliary facts 

4.1. Pricing options with digital and call options. In [IB], the author 
considered a certain market model where the bond is given by an increasing 
function Bt with Bq = 1, and we are given the underlying asset Xt and the 
equivalent martingale measure Q. The following class of payoff functions 
was considered. 

Definition 4.1. For a function f : IR + -4 IR, we denote f G n<g(Xr) i/f/ie 
following conditions are satisfied: 

(1) / is continuously differentiable except on at most countable set of 
points < si < S2 < • • • < Sat (and possibly on sq = 0) in which f 
and f have jump- discontinuities, 

(2) f(X T ) G L\Q), 

(3) / satisfies 

(4.1) lim \ f(x-)\Q(X T > x) = 
and, 

(4) t/^e Riemann-Stieltjes integral 

/•oo 

/ /'(a)y c (da) 

JO 

exists, where 

(4.2) F c (a) = E Q [ J B T 1 (X T -a) + ] 

denotes the price of European call option with strike a under measure 

Q. 

For Barrier type options, the class Hq(Xt,Y, A) of payoff function un- 
der consideration was a class of functions / such that the assumptions of 
Uq(Xt) was satisfied, where on condition (|4.ip . Q(Xt > x) was replaced 
with Q(Xt > x,Y G A) for some Borel measurable set A and V (a) was 
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replaced with 



V c > Y ' A (a) = m Q [B-\X T - a)+l YeA \. 

One of our main tools are the following pricing formulas. For proofs and 
details, we refer to [15]. 



N 

+ B T 1 Y,&-f(sk)Q(X T >s k ) 

k=0 
N 



Theorem 4.1. Let f E Hq(Xt)- Then the price of an option f(Xx) is 
given by 

/•oo 

V* = B^HO) - / f'(a)V c (da) 
Jo 

(4.3) 

+ B T 1 J2^+f(sk)Q(X T > s k ), 

k=0 

where A_/(sfc) = /(s&) — /(sfc— ) and A + /(s/ c ) is a jump on right side at 
Sk, respectively. The jump from the left at zero is defined as A_/(0) = 0. 

Corollary 4.1. Let f £ Hq(Xt)- If f is absolutely continuous on every 
interval (sfc, s fc+i)j then 

/■OO 

Vf = B^f(0) + / f"(a)V c (a)da 
Jo 

N 

+ B T 1 Y J ^-f(sk)Q(X T >s k ) 

k=0 
N 

+ B T 1 Y J ^+f{sk)Q(X T > s k ) 

k=0 

N 



(4.4) 



fc=0 

Corollary 4.2. If f is a linear combination of convex functions, then 

/•oo 

(4.5) Vf = B? l f(0) + /;(0)E Q (Xr] + / V c (a)n(da), 

Jo 

where [i is the measure associated with the second derivative of f . 
For barrier type options, the following theorem was proved in |15j . 
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Theorem 4.2. Let f £ Uq(Xt, Y, A). Then the price of the barrier option 
f(X T )l Y £A is given by 

POO 

V^ A = B~ l f(0)Q(Y £ A) — f'(a)V C ' Y ' A (da) 

Jo 

N 

(4. 6 ) +B T 1 J2&-f(sk)Q(X T >s k ,YeA) 

k=0 

N 

+ B T X & + f(s k )Q(X T >s k ,Ye A). 

k=0 

Corollary 4.3. Let f E Uq(Xt,Y, A) and assume that f'(a) is absolutely 
continuous on every interval (s k , Sfc+i)- Then the price of the barrier option 
f{X T )l Y aA is given by 

roc 

V f> Y > A = fl-V(O) + / f"(a)V C > Y ' A (a)da 
Jo 

N 

+ B T X &-f(s k )Q(X T >s k ,Ye A) 

k=0 

(4.7) N 

+ B T l Y &+f(s k )Q(X T >s k ,YG A) 

k=0 

N 

+ " f'{s k -))V c ^ A {s k ). 

k=0 

4.2. Rate of convergence for call and digital options. The following 
error expansions in a generalised binomial model for digital and call options 
was proved in [3]. 



Theorem 4.3. Let At = — and X„ 

n n 

the n-period binomial model, where 
u _ e aVM+\„cr 2 At 



an arbitrary bounded function of n. For 



-ay/At+\ n a 2 At 



if the initial stock price is Sq and the strike price is a and maturity T, then 
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(1) the price of a digital call option satisfies 



(4.8 



V^ (ft) (o) = V^(a) + 



+ o 



-rT - — ^ 

e e 2 



An _ d 2 (A n ) 2 | £ n ' 
2n n 



and 

(2) the price of a European call option satisfies 



(4.9) 



BT(n)V y 24(j V2^T n 



+ o 



1 



nWn 



where A n , A n and B n are given in Appendix A. 



Remark 4.1. In [3], the authors actually stated their result in a form where 
the remaining term was given by o(n _1 ). However, by careful examination 
of the proof one obtains that it is actually O ■ 

Remark 4.2. In the proof of Theorem \4-3\ the authors assumed that the 
strike a satisfies Sod n < a < ^Sou 71 . This is a natural assumption since in 
the n-period binomial model, the price of a digital option is e~ rT and the 
price of a call option is Sq if a < Sod n . Similarly, for a > Sou n the prices 
for both are zero. However, it is straightforward to see that the formulas of 
Theorem \4-3\ are true for any strike a. 

Remark 4.3. In [3j, the authors used payoff function f(x) = l x > a for 
digital option with strike a. Evidently, if the strike a is not on a lattice 
point, the formula |^75|) holds for the payoff f{x) = l x > a as well. For our 
result, we need similar formula for payoff f{x) = l x > a when the strike is on 
a lattice point. By examining the proof of Chang and Palmer |3J, we obtain 
that in this case, the formula is given by |^,<g| j ; where A n is replaced with 
A n (a) — 2. For the details, see [3]. 
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The following theorem for error expansion in CRR model for barrier call 
options was derived in |12j . 



Theorem 4.4. Suppose the initial stock price is So in the n-period CRR 
model. Then the asymptotic formulas for the prices of European barrier call 
options with strike a, barrier H , and maturity T are as follows: 



(1) for down-and-out (do) and down-and-in call (di) options with So > 
H: 

(a) if H < a, or H = a and H = a is not a terminal stock price: 

rC,do l N irC,do t - - ■ ■ IT, 1 



vS?„ ) («) = yS to w + [^]^ 

+ [b x - D X (A n ) 2 - Et (A^) 
1 



/A 2 



+ o 



■n 



V^ d i n) (a)=Vgf(a)-[A 1 A*] 



n 



l - + o 

n 



n x n 



1 



n\/n 



- B 2 + D 2 (A„) 2 - Ei (A^ 
(b) if H > a, or H = a and H = a is a terminal stock price: 

T C,do , \ -<rC,do, > r . . 1 



^T(°n)(«) = ^T(«) + [^]^ 



+ 



^3 - Ce 2 n - E 2 (A^) 2 



n 



n V nJn 



1 



+ o 



n 



1 



n\/n 



- B 4 + D 3 (A n ) z - Ce z n - E 2 (A£ 

(2) for up-and-in (ui) and up-and-out options (uo) with So < H, if 
H > K: 

1 



+ 



B 5 -D 2 (A n ) 2 -Ce 2 n -E 3 (A%y 



l - + o 

n 



1 



16 



vSj,w=ysr(«)-[^]^ 



r C.uo 
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- + o 

n \ ri\/n 



B G + Dx (A n ) 2 - Cel - E 3 (A* ) 



where &5 , A n , Aj, £>j, C, -Dj, j5j and e n are given in Appendix B. 



5. Proofs 



The proof of Theorem 12.11 is based on the following simple lemma. 



Lemma 5.1. Let f E II. Then 



(5.1) 

/or every p > 1 . 



POO 

/ |a/'(c 

JO 



e 2 (if 



da < oo 



Proof of Theorem \2.1\ Let n be fixed and assume first that / is continuous. 
We stress the fact that in both models, the derivative of the price of the call 
option with respect to the strike price is the price of negative digital option. 

Evidently, we may apply Theorem 14.11 for pricing option in n-period bi- 
nomial model. In order to apply Theorem 14.11 in Black-Scholes model, it is 
sufficient to show that 



(5.2) 



f'(a)Vg s (a)da 



< oo, 



where Vg S (a) = e rT $(d2). Note that for every number x > 0, we have 

1 

1 - <&(z) < ; e 2 . 



2-irx 



Hence for sufficiently large a, we have 

1 



(5.3) 



$(d 2 ) < 



Hence, by Lemma 15-H we have 15.21 
By ([4^8]) we have 



A n (a) I 

n n 



e~^~ < \di\ s e~^. 



Vgr {n) (a)-Vg 3 (a) 



B n (a) 



d 2 A 2 n (a 



+ 



n,\ n 



O n {a) 
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for some function O n {a) which is bounded in n. Thus, by Theorem 14.11 
obtain 



1 f°° 



m do 

— rT - 

e e 2 



1 f°° 
+ ~ / /'(«) 



2vr 

-rT 

e e 2 



-A„da 



+ 



1 



Note that for sufficiently large a, we have, 



2?r 

f'(a)O n (a)da. 



B r , 



da 



(5.4) 
and 
(5.5) 



e 2 



< ce 2 



d 2 d 2 

e~^~ |A n | < ce~~|di| 



for some constant c. Hence the integrals are finite and bounded in n. 
To complete the proof, we show that 

1 /*°° 

\f'(a)O n (a)\da = o(n- 1 ). 



Note that O* =: sup a sup n |O n (a)| < oo. We have 

1 



\f'(a)O n (a)\da 



1 



+ 



i 



\f'(a)O n (a)\da 
\f'(a)O n (a)\da 



< 



O* sup 0<a<1 | /'(a) | 



+ 



n^Jn 

oo 



\f'(a)O n (a)\da. 
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Evidently, the first term is o(n _1 ). For the second term, we obtain by Holder 
inequality that 
1 



< 



< cfe 



\f'(a)O n (a)\da 

roc 





1+e 


[ l+e d « 




J i a e 





1+t \O n (a)\ 1+e , ' 
all - — , , da 



n |(i+e) 



i 

l+e 



(o*y 

n< 



sup / | af (a) | 

n JO 



l+e 



| On (a) 



■da 



i 



where 7 = 2(e+i) > 1 if e > 2. It remains to show that 

a+e\°n(a)\ 



sup / I af' (a) I 

n JO 



-da < 00. 



ri\/n 



By ()5.ip . (|5.4p and (|5.5[) this is equivalent to condition 

\af'(a)\ 1 ^V£ T{n) (a)da< 



sup 
n Jo 



00. 



But this follows from the fact that binomial price converges to Black-Scholes 
price. Finally, by same analysis and Remark 14.31 we obtain that the sum- 
mation terms which correspond the jumps of / converges and the remaining 
term is o(n _1 ). This completes the proof. □ 



Proof of Theorem \2.2L The proof is essentially the same as the proof of The- 
orem [2J] and the details are left to the reader. □ 



Proof of Theorem \2.3l The proof relies on Theorem 14.41 Corollary 14.31 and 
the arguments in the proof of Theorem 12. 1[ We omit the details. □ 
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Appendix A. List of functions and constants for stock options 



u = e <rV~Ai+X n At g _ e -cry/Ai+\ n At 



A n = 1 — 2frac 



log ^ + n log d 



B„ — 



d\ + d x d\ + 2d 2 - 4di (2 - d x d 2 - dj)VT . 2 , 
H - i (r-A n fJ z ) 



24 

i \ 2\2 



6a 



A n = -a 2 T{6 + d\ + d 2 2 ) + AT{d\ - d 2 2 ){r - \ n a 2 ) - \2T 2 {r - X n a 2 ) 2 , 



Jn. — 



e 2 



2vr 



e 2 



2vr 



A ra d 2 A 2 | ff ra 
yfn 2n n 



2 _ d 2 (A ra - 2) 2 

2n n 



[A n -12a 2 T(A 2 -l)] 



24av / 2^T 

Appendix B. List of functions and constants for Barrier 

OPTIONS IN CRR MODEL 



u 



„oVAt 



-aVAt 



Ih = 



A 



H 



2frac 



log ^ + n log d 



2(7^ 2 

frac(2lff), for down options 
frac(—2ln), for up options. 

0, if the effective barrier is not a terminal stock price 

1, if the effective barrier is a terminal stock price. 
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dn ^° g (f) + (r + ^)r d ^ dn _ aVf 

o\JT 

d 21 = — — f= ; d 22 = d 21 - aVT 

o\JT 

log (f ) + (r + lg 2 ) T 
"31 = /= 5 "32 = "31 - ffVi 

cry J 

lo g(f) + ( r + ^ 2 ) T 

C?41 = p= 5 d,42 = ^41 - (xVT 

1 2 

r - ^o A A cr 

cr 4 - 4<r 2 r + 12r 2 - _ err 

9i = 2T (a 2 d l2 + /?Vt) + _ ^ (1 _ 4), i = 1, 2, 3, 4 

5, = 2T (a 2 ^ + /?Vt) + - ^ (1 - 4), i = 1,2,3,4 

Gi = ^e ibn-m), G > = ^{h) e '(92 - 92) 

Ai = 4Vf/»l(d2l,d22), ^2 = 4Vr/li(d41,d42) + ^ e "^ " f ) 

a 3 = 4Vf /M(-d4i, -d 42 ) - 4Vr/n(-d2i, -da 2 ) - ^ e "^ (i - ^) 



where 



2r 

M*,y)=(§) ^ (h&^x) - ^LHa^(y) J , i = 0,1,2. 

B 1 = G 1 -G 2 + Ih (d 21 ,d 22 ), B 2 = B 1 -G 1 
B3 = G3 — G4 + Iho(d4i, c/42), B4 = B% — Gi 

-B5 = G2 + G3 — G4 + Iho(—d 2 i, —d 22 ) — Iho(—d4i, —^42), B$ = B§ — G\, 
where 

/4/3 + fa 3 V__7<So 



log ( ^ ) r. 



Ci 

c 3 

Di 

E3 



F, 



n,l 



n,2 



n,5 



n,6 



n,7 
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c 2 — C3 
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UJ C2 = 7^ e ( d31 "^ 32 ) 

s 



S _is_L (, a \ 



. 2tt 

8TH 2 (d 21 ,d 22 ) + Ci, £ 2 = 8Th 2 (d 41 ,d 42 ) + 3C>2 2 +C3 
8Th 2 (-d 21 , -d 22 ) - 8Th 2 {-d 41 , -d 42 ) - d + 3C2 + C \ 



A = -a 2 T(6 + d\ + d\) + 4Tr{d\ - d 2 2 ) - 12r 2 T 2 , 



_ *L 

S e 2 



[A-12a 2 T(A 2 -l)], 



n,3 



-^1,4 = 



.4 2 A^, a < # 
AiA%, a>H 

B 3 -Ce 2 n -E 2 (A«) 2 , a<H 
B 1 -D 1 (A n ) 2 -E 1 (A») 2 , a>H 

-B 4 - D z {A n ) 2 + Ce 2 n + E 2 {A*)\ a<H 
-B 2 - D 2 (A„) 2 - E 1 {A*) 2 , a>H 

A 3 A%, a<H 

0, a>H 

B 5 -D 2 (A n ) 2 -Cel-E 3 (A»y 
H, 

-B 6 -D 1 (A n ) 2 + Ce 2 n + E 3 (A%) 2 , a<H 

1, a>H 

-B 6 -D 1 (A n ) 2 + Ce 2 n + E 3 (A^) 2 , a<H 
0, a>H 



a<H 
a> H 
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